Truncated matricial moment problems on a finite 
interval: the operator approach. 



S.M. Zagorodnyuk 



1 Introduction. 



In this paper we study the following problem: to find a non-decreasing 
matrix function M{x) = (W'fc,«(^))fc^/=o [^j^]) which is left-continuous in 
{a,b),M{a) 



0, such that 



x^dM{x) = Sn 



n 



0,1, ...J, 



(1) 



where {Sn}n=Q is a given sequence of Hermitian {N x N) complex matrices, 
iV G N, / e Z+. Here o, 6 G M: a < b. 

In the scalar case this problem was solved by M.G. Krein, see p]. Re- 
cently, a deep investigation of the matrix moment problem ([T|) was completed 
by A.E. Choque Rivero, Yu.M. Dyukarev, B. Fritzsche and B. Kirstein, 
see These authors used the Potapov method for interpolating prob- 

lems which was enriched by the Sachnovich method of operator identities. 
Set 



( So 
Si 



i+j)i,j=0 



Si 

S2 



\ Sk s , 



k+l 



Pfc = {-abSi+j + (a + b)S, 



i+j+i 



S 



Sk \ 

Sk+i 

S2k / 



\k-l 



k e 



k G 



2k < I; (2) 



2k < I. (3) 



If we choose an arbitrary element / = (/o, /i, . . . , Jn-i), where all fk are 
some polynomials and calculate J^^ fdM f* , one can easily deduce that 

Tk>0, k eZ+ : 2k <l. (4) 

In the case of an odd number of prescribed moments / = 2d, the strong result 
of A.E. Choque Rivero, Yu.M. Dyukarev, B. Fritzsche and B. Kirstein is that 
conditions 

Td >0,fd> 0, (5) 

are necessary and sufficient for the solvability of the matrix moment prob- 
lem ([1]), see jSl Theorem 1.3, p. 106]. For the case > 0, P^ > 0, they 
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parameterized all solutions of the moment problem via a linear fractional 
transformation where the set of parameters consisted of some distinguished 
pairs of meromorphic matrix- valued functions. 
Set 

Hk = {—aSi^j+Si^j^i)'^ j^Q, Hk = {bSi^j — Si^j^i)'lj^Q, G Z+ : 2k+l < I. 

(6) 

In the case / = 2d + 1, the analogous to the result for / = 2d, the result of 
A.E. Choque Rivero, Yu.M. Dyukarev, B. Fritzsche and B. Kirstein states 
that conditions 

Hd > 0, Hd> 0, (7) 

are necessary and sufficient for the solvability of the matrix moment prob- 
lem (P, see [21 Theorem 1.3, p. 127]. For the case Hd > 0,Hd > 0, they 
parameterized all solutions of the moment problem via a linear fractional 
transformation. The set of parameters consisted of some distinguished pairs 
of meromorphic matrix- valued functions. 

In this work we will study the matrix moment problem ([1]) by virtue of 
the operator approach based on the use of the generalized resolvents of 
some symmetric operators. In the study of the classical Hamburger moment 
problem this approach finds its origin in the papers of M.A. Neumark 
and M.G. Krein, M.A. Krasnoselskiy [6], see also [7]. All these authors used 
orthogonal polynomials connected with a Jacobi matrix related to the mo- 
ment problem. Lately, we showed that this method in a general setting can 
be applied to the Hamburger moment problem both in the non-degenerate 
and degenerate cases, see [8]. Our goal here is to describe all solutions of 
the matrix moment problem ([T]) in a general case. This means that no con- 
ditions besides solvability of the moment problem will be assumed. At first, 
we study the case of an odd number of prescribed moments / = 2d, and 
then we shall reduce the case of an even number of moments / = 2d + 1 
to the previous case {d € In our study we shall use the basic results 

of M.G. Krein and I.E. Ovcharenko on generalized sc-resolvents of symmet- 
ric contractions, as well as M.G. Krein's theory of self-adjoint extensions of 
semi-bounded symmetric operators, see |9].|10j.[TT|. 

Notations. As usual, we denote by M, C, N, Z, Z+ the sets of real, complex, 
positive integer, integer, non-negative integer numbers, respectively. The 
space of n-dimensional complex vectors a = (ao, ai, . . . , a^-i), will be de- 
noted by C", n G N. If a € then a* means the complex conjugate vector. 
For a complex (n x n) matrix A, we denote by Ker A a set {x € C" : Ax = 0}. 
By P we denote a set of all complex polynomials and by ¥d we mean all com- 
plex polynomials with degrees less or equal to d, d G (including the zero 
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polynomial). Let M{x) be a left-continuous non-decreasing matrix function 
M{x) = Q on R, M{-oo) = 0, and tm{x) := Y.k=o '^k,k{x); 

^'(x) = {dm}^^i/dTM)khlQ- We denote by Lp'{M) a set (of classes of equiv- 
alence) of vector functions / : M — > C^, / = (/o,/i, • • • ,fN-i), such that 
(see, e.g., [T2]) 

WfWhiM) ■■= I f{x)^{x)nx)dTM{x) < oo. 

The space Lp'{M) is a Hilbert space with the scalar product 

(/,5)l2(m) := / f{x)^{x)g*{x)dTM{x), f,g G L\M). 

For a separable Hilbert space -ff we denote by {■,-)h and || • \\h the scalar 
product and the norm in H, respectively. The indices may be omitted in 
obvious cases. 

For a linear operator ^ in if we denote by D{A) its domain, by R{A) its 
range, and by A* we denote its adjoint if it exists. If A is bounded, then 
\\A\\ stands for its operator norm. For a set of elements {xn}neB in H, 
we denote by Lm{xn}neB and span{2;„}„gB the linear span and the closed 
linear span (in the norm of H), respectively, where B is an arbitrary set of 
indices. For a set M Q H we denote by M the closure of M with respect to 
the norm of H. By Eh we denote the identity operator in H, i.e. Ehx = x, 
X £ H. If Hi is a subspace of H, by Ph^ = we denote the operator of 
the orthogonal projection on Hi in H. A set of linear bounded operators 
which map H into H we denote by [H]. 

2 The case of an odd number of given moments: 
solvability and a description of solutions. 

We shall use the following important fact (see, e.g., [13^ p.215]): 

Theorem 2.1 Let K = (Kn,m)nm=o ^ 6e a positive semi-definite com- 
plex ((r-|-l) X (r-|-l)) matrix, r £ Z_|_. Then there exist a finite- dimensional 
Hilbert space H with a scalar product (•, •) and a sequence {xnYn=o ^ ' 
such that 

Kn,m = {Xn,Xm), 71, m = 0, 1, r, (8) 

and span{x„}^^o = ^■ 
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Proof. Let {a;„}^^Q be an arbitrary orthonormal basis in C". Introduce 
the following functional: 

r 

[x,y\= ^ Kn,manbm, (9) 
n,m=0 

for x,y G C", 



r 

X = 

n=0 m=0 



The space equipped with [•, •] will be a quasi-Hilbert space. Factorizing 
and making the completion we obtain the required space H (see, e.g., [14:\ 
p. 10-11]). □ 

Consider the matrix moment problem ([1]) with / = 2d, d E N. Suppose 
that > (as we noticed in the Introduction, condition Q is necessary 
for the solvability of the moment problem). Let Td = {7d;n,m)^^lo ^, 
ld;n,m ^ C. By Theorem 12.11 there exist a finite-dimensional Hilbert space 
H and a sequence {xn\^n=o^^ ^ ™ such that 

{Xn,Xm) = 'ld;n,Tn, n, m = 0, 1, (d l)iV - 1, (10) 

and span{x„}^'!^o^^^ ^ = Linjxnj^'ty^^ ^ = H. Notice that 

7d;rN+j,tN+n = Sr+f < j, 71 < iV - 1; < r,t < d, (11) 

where 

are the given moments. From (|lip it follows that 



ld;a+N,b = ld;a,b+N, « = rN+j, b = tN+n, <j,n < N-1; < r,t < d-1. 

(12) 

In fact, we can write 

7d;a+N,b = 7d;{r+l)N+j,tN+n = ^r+t+l = ld;rN+j,{t+l)N+n = ld;a,b+N- 

Set Ha = {a^n}^=o ^- We introduce the following operator: 

dN-l dN-1 
Ax= UkXk+N, X € Ha, X = ^ OfcXfc. (13) 

A:=0 A:=0 

The following proposition shows when the operator A is correctly defined. 
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Theorem 2.2 Let a matrix moment problem (1^ with I = 2d, (i € N, be 
given. The moment problem has a solution if and only if conditions (0j are 
true and 

KerTd-i CKerfrf.i, (14) 

where Td-i = {Si+j+2)'!;jto- 

If conditions [Tl\ ) are satisfied then the operator A in [T3\) is correctly 
defined and the following operator: 

Bx = -^A - ^Eh, X G Ha, (15) 
— a — a 

is a contraction in H (i.e. \\B\\ < 1). Moreover, operators A and B are 
Hermitian. 

Proof. Let the matrix moment problem ([T]) has a solution M(x) = 
i'mk^i{x))^f2o- Consider the space Lp'{M) and let Q be the operator of 
multiplication by an independent variable in Lp'{M). The operator Q is 
self-adjoint and its resolution of unity is (see [12j ) 

Eb-Ea = E{[a, b)) : h{x) ^ X[a,b){x)h{x) , (16) 

where X[a,b){x) is the characteristic function of an interval [o, 6), — oo < a < 
b < +00. 

Set 4 = {ekfi,ek,i,. ■ ■ ,ek,N-i), e^j = < j < iV - 1, for A; = 

0, 1, . . . A'' — 1. A set of (classes of equivalence of) functions / E Lp'^M) 
such that (the corresponding class includes) / = (/o, /i, . . . , /at-i), fj S P^, 
we denote by P^(M) and call a set of vector polynomials of order d in 
L2(M). Set i^o(^) = Since P^(M) is finite-dimensional, we have 

Llo{M)=¥l{M). 

For an arbitrary polynomial (in a class) from P^(M) there exists a unique 
representation of the following form: 

N-l d 

fix) = ^ '^akjx^ik, Okj € C. (17) 
k=o j=0 

Let a polynomial g G P^(M) have a representation 

N-l d 

9ix) = T.T.f^^'rx'ei, Pi,reC. (18) 

1=0 r=0 
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We can write 

N-l d „ N-1 d 

{f,9)L2{M) = X] X] "fcjA:^ / x^^''ekdM{x)el = X] X] ^k,jl3i,r* 

k,l=Oj,r=0 ''^ k,l=Oj,r=0 

„ N-l d 

* / x^^''dmk,i{x) = X] X] '^k,jPi,rsfl^, (19) 



A:,/=0 j,r=0 



where 



•^n = i4Xr=o^ " ^ (20) 

are the given moments. On the other hand, we can write 

d N-l d N-l \ N-l d 

CtkjXjN+k,''^ ^ Pl,rXrN+l 1 = ^ ^ 0!kjPl,r{XjN+k, XrN+l)H 
^j=0 k=0 r=0 1=0 J ^ k,l=Oj,r=0 

N-l d N-l d 

= X/ X/ ^kjPl,rld;jN+k,rN+l = X/ X/ (^1) 
fc,/=Oj,r=0 k,l=Oj,r=0 

From relations ()19p . (|2ip it follows that 

(d N-l d N-l \ 

X] ^ 0!k,jXjN+k, ^ ^ (3l,rXrN+l 1 • (22) 
j=0 fc=0 r-=0 1=0 J ^ 

Set 

d N-l 

^/ = X] X] ak,jXjN+k, (23) 

j=0 fc=0 

for /(x) E P2(^)^ j(^) ^ E£o' E'=o«fcJ^'efc, «A:,i G C. 

If /, g have representations (fT7l) . (fT8| . and ||/ — fl'||L2(^) = 0, then from ([22]) 

it follows that 

\\Vf-Vg\\l = (Vif-g),V{f-g))H = {f-9,f-9)LHM) = \\f-9\\h^M) = 0- 

Thus, y is a correctly defined operator from P^(M ) to H. 

Relation (I22p shows that V is an isometric transformation from P^(M) onto 

hm{xn}fj^Q~^^ ^ ■ Thus, V is an isometric transformation from q{M) onto 
H. In particular, we note that 

Vx^ik = XjN+k, 0<j<d; 0<k<N -1. (24) 
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Set Lj iiM) := L'^{M) Q LI q{M), and U := V ® ^(jv/)- The operator U 
is an isometric transformation from Lp'{M) onto H ® L'^ =• H. Set 

A:=UQU-^. (25) 
The operator ^4 is a self-adjoint operator in H. Notice that 
UQU~^XjN+k = VQV'^XjN^k = VQx^Sk = Vx^+^Sk = xq+i)^^^ = Xj^+k+N, 
0<j<d-l; < k < N - 1. 

By hnearity we get 

dN-l dN-1 

UQU'^x = ^ akXk+N, X e Ha, X = OkXk- 

k=0 k=0 

Consequently, the operator A in (fT3|) is correctly defined and 

A = A\h^. (26) 
Since A is correctly defined, from the equality 

^kXk = 0, (27) 

fc=0 

with some complex numbers ^fc, it should follow the equality 

dN-l 

^kXk+N = 0. (28) 

k=0 

On the other hand, the equality ([27|) is equivalent to the equalities 

dN~l dN-l 

Y ^k{xk, Xl) = Y ^kld;k,l = 0, / = 0, 1, ...,dN - 1. (29) 
k=0 k=0 

Analogously, the equality ([28|) is equivalent to the equalities 

dN~l dN~l 

Ck{xk+N,Xl+N) = Y ^kld;k+N,l+N = 0, / = 0, 1, dN - 1. (30) 
fc=0 A,-0 

If we shall use the matrix notations, the equality 

\dN~l 



{Co,Cl,---,^dN~l){ld;k,l)kl=0 =0' (31) 



implies the equality 

iS,0, (.1, (,dN-l){7d;k+N,l+N)i^=0 = 0. (32) 



Thus, relation (jl4p is true. 
Consider the following operators: 



R ■■- 7 Q - t-^—El2,m), (33) 

— a — a ^ ' 

B := URU-^ = -^A - ^En. (34) 
b-a b- a ^ 

Define an operator B by the equality (fT5|) . From ([26l) . ([M]) we get 

B = B\h^. (35) 
For an arbitrary / E D{R) = D{Q) we can write 



\\Rff 



2 a + b 

-X 



b — a b — a 



f{x)dM{x)r{x) < / f{x)dM{x)r{a 
J a 



= \\ff, 

and therefore the operators R, B and B are contractions. Since R is Her- 
mitian, the operators B, B are Hermitian, as well. Choose an arbitrary 
x G Ha, X = Ylk=o^ otkXk, and write 

dN~l dN-l 

< \\x\\'^ — \\Bx\\'^ = akaj{xk,Xj) — akOjiBx^, Bxj) = 

k,j=0 k,j=0 
dN-l dN-l / „ , i, n , u \ 

X - / 2 a + b 2 a + b \ 

akaj^d;k,j- > , OfcOj 7 Xk+N - 7 Xk, Xj+N - T Xj = 

^ \b- a b-a b-a b - a J 

—( 4 2(a + 6) 2(a + 6) 

= 2^ afcOj \^d:k,j - _ ^y 7d;k+N,j+N + "^^^^Irfifc+iVj + J^^^^ld;k,j+N- 

-Jb^^'-'^'O ■ ^''^ 

If we multiply the both sides of the latter inequality by (6 — a)^ and use (fT2|) 
we get 

dN-l 

< ^ akOj ((6 - af^d;k,j - A'yd;k+N,j+N + 4(a + b)'Jd;k+N,j - (a + bf-fd;k,j) = 
k,j=0 
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dN-l 

= ^ akOji-Aab'jd-kj - ^ld;k+N,j+N + 4(a + b)-^d;k+N,j)- 

k,j=0 

Therefore < {-ab-fd-k,j - 7d;k+N,j+N + (a + b)-fd-k+N,j)tjZo = ^rf' ^nd 
the second relation in ([5]) is true. The necessity of the first relation in ([5]) 
was discussed in the Introduction. Thus, we established the necessity of 
conditions ([5l). (jl4p for the solvability of the moment problem ([T|). 

Let a matrix moment problem ([T]) be given and conditions (l5]). (ll4p be 
true. We construct a Hilbert space H and a sequence {xn}^=Q^^ ^ in H, 
such that relation pOj) is true and span{xn}^'^o^^'^ ^ = H. Condition (fT^ 
provides that equality (j27p will imply equality (128p . This means that the 
operator A in (jl3p is correctly defined. We define an operator B by the 
equality ([H]). For arbitrary x,y £ Ha = span{x„}^^o\ x = Ylk=o^ '^kXk, 
y = Yl'j=o^ Pj^ji using ([12]) we can write 

dN-1 dN-1 

{Ax,y) = ^ ak(3j{xk+N,Xj) = ^ akPj{xk,Xj+N) = {x, Ay). 

k,j=0 k,j=0 

Thus, operators A and B are Hermitian. If we use relation (j36p (except the 
first inequality in it) and the second condition in ([5]), we obtain that B is 
a contraction. By Krein's theorem |1H Theorem 2, p. 440], there exists a 
self-adjoint extension B of the operator B in H with the same norm as B 
(and therefore it is a contraction). Let 

B = j'^XdE^, (37) 

where {-Ea} be the left-continuous in [—1, 1), right-continuous at the point 
1, constant outside [—1, 1], orthogonal resolution of unity of B. Choose an 
arbitrary a, < a < d{N a = rN + j,0<r<d, 0<j<N-l. 

Notice that 

Xa = XrN+j = ^X(^j._i-^]\fj^j = ... = A^Xj. 

Then choose an arbitrary /?, < /3 < d{N + 1) - 1, P = tN + n, < t < d, 
< n < N — 1. Using (jlip we can write 

S^+t = 1d;rN+j,tN+n = {XrN+j , XtN+n) H = {A'' Xj , A^ Xn) H = 

-^B + ^E^^ x„[^-^B + ^E^^ xJ = 
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h-a~ a + b^Y^' \ fWb-a, a + bV^'- 

-BH ^^hJ Xj,Xnj = J d{ExXj,x 

Set 

fhj^nix) = [E 2 ^ a+bXj,Xr,)H, < j,n < N - 1. 

b—a b—a 

Then 

rb 

sl'^^ = / x''+^dmj^n{x), 0<j,n<N-l, 0<r,t<d. (38) 

J a 

Prom relation (j38p we derive that the matrix function M(A) = ifnj^n{x))^~J;Q 
is a solution of the matrix Hamburger moment problem ([1]) (Properties of 
the orthogonal resolution of unity provide that M(A) is left-continuous in 
{a,b), non-decreasing and M(a) = 0). 

It remains to prove the last statement of the Theorem. If conditions ([5]). p4p 
are satisfied then we proved that the moment problem ([1]) has a solution. In 
this case we showed that the operator A in (|13p is correctly defined and the 
operator B in (jl5p is a Hermitian contraction in H. The fact that operators 
A and B are Hermitian was established, as well. □ 

We shall continue our considerations before the statement of Theorem l2.2[ 
We assume that conditions ([5|). (jl4p are true. Therefore the operators A 
in (|13p and B in psp are correctly defined Hermitian operators and \\B\\ < 1. 
Let B be an arbitrary self-adjoint extension of i? in a Hilbert space H ^ H. 
Let Rz{B) be the resolvent of B and {i?A}AeM be an orthogonal resolution 
of unity of B. Recall that the operator-valued function = Rz{B) is 

called a generalized resolvent of B, z ^ C\M. The function = P^Ex, 
A € M, is a spectral function of a symmetric operator B (e.g. [15J). There 
exists a one-to-one correspondence between generalized resolvents and (left- 
continuous or normalized in another way) spectral functions established by 
the following relation (|16j): 

{-RJ,g)H= [ Y^d{Bxf,g)H, f,g€H,zGC\R. (39) 

To obtain the spectral function from relation (j39p . one should use the 
Stieltjes- Perron inversion formula (e.g. [7j). 

In the case when i? is a self-adjoint contraction, the generalized resolvent 
R-^ = PhRz{B) is called a generalized sc-resolvent of -B, see [9], [10]. The 
corresponding spectral function of B we shall call a sc-spectral function of 
B. 
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Let ^ 

B = J ^XdEx, (40) 

where {Ex} be the left-continuous in [—1, 1), right-continuous at the point 
1, constant outside [—1,1], orthogonal resolution of unity of -B. In a similar 
manner as after (j37|) we obtain that 

fb 

sl'l^ = I x''+*(i%,„(x), < j, n < iV - 1, < r, t < d, (41) 

J a 

where 

mj^nix) = iPHE^^_a+bXj,Xn)H, 0<j,n<N-l. (42) 

b—a b—a 

Thus, the function M{x) = (f^j,n(2;))^^Q is a solution of the matrix moment 
problem ([1]). 

Theorem 2.3 Let a matrix moment problem ^^ with I = 2d, (i G N, be 
given. Suppose that conditions ^14\) are true. All solutions of the moment 
problem have the following form 

M{x) = {mj^n{x))f~}Q, mj.n{x) = (E^^._a+6Xj,Xn)H, < j,n < iV-1, 

h — a b — a 

(43) 

where is a left- continuous in [—1,1), right- continuous at the point 1, 
constant outside [—1, 1] sc- spectral function of the operator B defined in i fi5|) . 
Moreover, the correspondence between all solutions of the moment problem 
and left- continuous in [—1, 1), right- continuous at the point 1, constant out- 
side [—1,1] sc- spectral functions of B in JT^j ) is one-to-one. 

Proof. Choose an arbitrary left-continuous in [—1,1), right-continuous at 
the point 1, constant outside [—1, 1] sc-spectral function of the operator 
B from (jl5p . This function corresponds to a left-continuous in [—1, 1), right- 
continuous at the point 1, constant outside [—1, 1] resolution of unity {E\} 
of a self-adjoint contraction B ^ B m. a Hilbert space H ^ H. Considera- 
tions before the statement of the Theorem show that formula ()43p defines a 
solution of the moment problem. 

On the other hand, let M{x) = {fnk^i{x))^Y2o be an arbitrary solution of 
the matrix moment problem ([1]). Proceeding like at the beginning of the 
Proof of Theorem 12.21 we shall construct a self-adjoint contraction B ^ B 
in a space H ^ H. Repeating arguments before the statement of the last 
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theorem, we obtain that the function M{x) = {mj^n{x))^n=o^ where mj^nix) 
are given by (j42p . is a solution of the moment problem. 
Choose an arbitrary z G C\M and write 



j -^^—^d{ExXk,Xj)g=(^j YZ^dExXkjXj^ ^ = 
C -^dExXk, U-'x,) =( C ^dU-'ExUik, ej) 

[ j^dER,xek, e^-) = {{R - z^L2(A/))~^efc, ej) ^(m) = 

d{Euek,ej)L2^M) 



^ / a + b 



b — a b — a 



/■I 1 

= / Y^^*^-^— A+^i±^^'='^i)^''W' 0<A:,j<A^-l, (44) 

where {Ex\ and {E^-x} are left-continuous in [—1, 1), right-continuous at the 
point 1, constant outside [—1,1], orthogonal resolutions of unity of operators 
B = URU^^ and R, respectively. Here {Ex} is the orthogonal resolution of 
unity of Q, given by (jl6p . By the Stieltjes-Perron inversion formula we get 



{ExXk,Xj)g = {Eb^^^a±bek,ej)L2(^M), <k,j < N -I, (45) 

for each A G [— 1, 1] , such that A is a point of continuity for Ex and Eb^ \,2d± ■ 
Using the change of variable we obtain that 



frikjix) = (Ph E^^ a+bXk,Xj)H = (-E:c4,ei)L2(M)) (46) 

b—a b — a 



where < A:, j < A^ — 1, for x E [a, b]: x is a point of continuity of E 2 _a±b 

b—a b — a 



and Ex- Using ()16p we can write 

{Exek,ej)L2(^M) = rukjix), x G (a, 6), 

and therefore 



fhk,j{x) = mkj{x), (47) 
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where x G [a, x is a point of continuity of £'_2_ ._a±b and E^. Since 

b—a b — a 

matrix functions M[x) and M{x) are left-continuous in (a, 6), they coincide 
in (a, h). It remains to note that M{a) = M{a) = 0, and M{h) = M{h) = Sq, 
to obtain 

M{x) = M{x), X G [a,b]. 

Consequently, all solutions of the truncated moment problem are generated 
by left-continuous in [—1, 1), right-continuous at the point 1, constant out- 
side [—1,1] sc-spectral functions of B. 

It remains to prove that different sc-spectral functions of the operator B 
produce different solutions of the moment problem ([1]). Suppose to the 
contrary that two different left-continuous in [—1, 1), right-continuous at 
the point 1, constant outside [—1,1] sc-spectral functions produce the same 
solution of the moment problem. That means that there exist two self- 
adjoint contractions Bj ^ B, in Hilbert spaces Hj 5 H, such that 

P^^i?l,A/P^^i?2,A, (48) 

{P^^Ei^XXk,Xj)H = {PH'E2,xXk,Xj)H, 0<k,j<N-l, AG [-1,1], 

(49) 

where {En^\}xeR are left-continuous in [—1, 1), right-continuous at the point 
1, constant outside [—1, 1], orthogonal resolutions of unity of operators B^, 
n = 1,2. Set Ltv := Lm{xk}k=o,N-i- By linearity we get 

{Plj'E^^xx,y)H = iPS'E2,xx,y)H, x,yeLN, AG [-1,1]. (50) 

Denote by Rn,x the resolvent of Bn, and set R^.a '■= PH'^Rn,x, = 1)2. 
From ([50]), ([39]) it follows that 

{'Ri,xx,y)H = {'R2^xx,y)H, x,y£LN, A G C\M. (51) 

Choose an arbitrary z G C\M and consider the space := {B — zEH)Ha- 
Since 

Rj^,{B - zEh)x = {Bj - zEH,r^ {Bj - zEh, )x = x, x e Ha = D{B), 
we get 

Ri,zU = R2,zU e H, n G H^; (52) 
Ri,2?x = R2,zW, uGH^, ze C\M. (53) 

We can write 

(Rn,zX,u)H = {Rn,zX,u)H„ = {x,Rn,zU)Hn = {x,'Rn,zU)H, X G H, U G H-, 
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n = l,2, (54) 

and therefore we get 

(Ri,2X,u)h = (R2,z2;,n)//, x£H,ueHY. (55) 

Choose an arbitrary u G H, u = Ylk=o^^^ '^k^k^ G C. Consider the 
following system of linear equations: 

^±^ + z]dk = Ck, k = 0,l,...,N -1; (56) 
b — a J 

^ -dk-N-(P^ + z]dk = Ck, k = N,N + l,...,dN + N-l; (57) 



b — a \b — a 

where {dk}'l^Q^^^ are unknown complex numbers, z £ C\M is a fixed pa- 
rameter, a,b are from ([T]). Set 



4 = 0, k = dN,dN +l,...,dN + N -1; 

a + b 
b — a 



dk-N = ((i-^ + z] dk + Ck] , k = dN+N-l,dN+N-2,...,N; 



(58) 

For such numbers ^, all equations in (f57|) are satisfied. Equa- 

tions (i56]l are not necessarily satisfied. Set v = Yl'l=o^~'^ dkXk = YUkl^o^ dkXk- 
Notice that v G Ha = D{B). We can write 

{B - zEh)v = (t^A - ^:^Eh - zEh] V = 
\b — a b — a ) 

f 2 fa + b \ 



A:=0 

dN+N~l 



b — a \b — a 



= g [b^/'-^-\b^a^')V'''^ 
where d_i = d-2 = ... = d-N = 0. By the construction of d^ we have 

{B - zEh)v - = ^- {^~^ + z^dk- Ck}j Xk] 

u = {B - zEh)v + X] ( ( I — ~ + z] dk + Ckj Xk, ue H, z e C\M. 

(59) 
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By ()59p an arbitrary element y H can be represented as y = yz + y', 
y- £ H-, y' G Ln- Using ([5T]) and (f55]) we get 

y)// = (Ri,^x,y2+y')H = (^2,^2;, y-+2/')H = {^2,zX,y)H, x G Lat, y G 
Thus, we obtain 

Ri,^x = R2,zx, X G Lat, z g C\M. (60) 

For an arbitrary x G L, x = Xz + x' , Xz G -ff^, G L^r, using rela- 
tions ([53]) ■ ([60|) we obtain 

Hl^zX = ^l^z{Xz + x) =Y{,2,z{Xz + x) =Ii2,zX, X G L, z G C\M, (61) 

and 

Ri,^x = R2,^x, X G i?, z G C\M. (62) 

By (I39p that means that the corresponding sc-spectral functions coincide 
and we obtain a contradiction. □ 

We shall recall some known facts about sc-resolvents (see [9], [10]). Let 
j4 be a Hermitian contraction in a Hilbert space H with a non-dense closed 
domain P = D{A)^ and TZ = H QT>. A set of all self-adjoint extensions of A 
in which are contractions, we denote by Bh{A)- A set of all self-adjoint 
extensions of A in a Hilbert space H ^ H, which are contractions, we denote 
by B^{A). The set BniA) is non-empty. Moreover, there are a "minimal" 
element and a "maximal" element A^^ in this set, such that Bh{A^) 
coincides with the operator segment 

A^'<A<A^^. (63) 

In the case A^^ = A^ the set Bh{A) consists of a unique element. This case 
is called determinate. 

In the case A^^ ^ A^'^ the set Bh{A) consists of an infinite number of ele- 
ments. This case is called indeterminate. 

The case A^^x 7^ A^x, x G 7^\{0}, is called completely indeterminate. The 
indeterminate case can be always reduced to the completely indeterminate. 
If 7^o = {x G 7^ : Ai'x = A^^x}, we can set 

AeX = Ax, x eV; AeX = A^'x, x G TZq. (64) 

The sets of generalized sc-resolvents for A and for A^, coincide ([10., P- 1039]). 
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Elements of Bh{A) are canonical (i.e. inside H) extensions of A and their 
resolvents are called canonical sc-resolvents of A. On the other hand, ele- 
ments of B^{A) for H D H generate generalized sc-resolvents of A. The set 
of all generalized sc-resolvents we denote by TZ'^{A). 
Set 

C = A^ - A^, (65) 
Q^{z)= (C^R'^C^ +Eh)\^, zeC\[-l,l], (66) 
where R'^ = {A^' - zEhY^. 

An operator- valued function k{z) with values in [H] belongs to the class 
Rn[-l,l] if 

1) k{z) is analytic in z G C\[— 1, 1] and 

\mkiz) ^ 
— — ^<0, zGC:Imz/0; 
Im z 

2) For z G ]R\[— 1, 1], k{z) is a self-adjoint positive contraction. 
Theorem 2.4 ( JT^ p. 1053]). The following equality: 

Rl = R'i- R^.C^ (Eh + (Q^(z) - E)k{z)r^ c'^Ri^, (67) 

where k{z) G i?7^[— 1,1], i?^ ^ '^'^(^); establishes a one-to-one correspon- 
dence between the set Rti[—1, 1] and the set TZ'^{A). 

Moreover, the canonical resolvents correspond in (67) to the constant 
functions k{z) = K, K ^ [0,En]. 

Comparing Theorem 12.31 and Theorem 12.41 we obtain the following result. 

Theorem 2.5 Let a matrix moment problem (1^ with I = 2d, d G N, be 
given. Suppose that conditions ^14^ are true. Let the operator B be 
defined by M5\) . The following statements are true: 

1) If = B^^ , then the moment problem has a unique solution. This 
solution is given by 



M{x) = (rnj,„(x))^„io, mj,n{^) = {E^J^^_a+bXj,Xn)H, 0<j,n< N-1, 

b — a b—a 

(68) 

where {Ez} is the left- continuous in [—1, 1), right- continuous at the point 1, 
constant outside [—1,1], orthogonal resolution of unity of the operator A^. 
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2) If 7^ B , we define the extended operator B^. (see the construction 
in ). All solutions of the moment problem have the following form 

M{x) = (mj,„(x)):^~io, "^j>(2;) = {^^^_a+bXj,Xn)H-, <j,n< N-1, 

b—a b — a 

(69) 

where is a left- continuous in [—1,1), right- continuous at the point 1, 
constant outside [—1,1] functions given by 

iR^,-R^,C-2 (Eh + iQ^^iz) - E)k{z)y' c"2R^J,g)H = [ -^diExf,g)H, 

f,geH, zeC\M, (70) 

where k{z) G i?7^[— 1,1]. Here IZ = H Q D{B(.), and C,Q^{z),Rz are de- 
fined by \65\ ). [6 ^) with A = Bg. Moreover, the correspondence between all 
solutions of the moment problem and k{z) € i?7^[— 1, 1] is one-to-one. 

Proof. It remains to consider the case 1). In this case aU self-adjoint 
contractions i? ^ in a Hilbert space H ^ H coincide on H with B^^, 
see [ini P- 1039]. Thus, the corresponding sc-spectral functions are spectral 
functions of the self-adjoint operator B^. However, a self-adjoint operator 
has a unique (normalized) spectral function. Thus, a set of sc-spectral func- 
tions of B consists of a unique element. This element is a spectral function 
o^B^'. a 

Consider a matrix moment problem ([T]) with 1 = 0. In this case the 
necessary and sufficient condition of solvability is 

So > 0. (71) 

The necessity is obvious. On the other hand, if relation (j7ip is true, we can 
choose 

M{x) = ^So, xG[a,b]. (72) 
b — a 

This function is a solution of the moment problem. A set of all solutions 
consists of non-decreasing matrix functions M(x) on [a,b], left-continuous 
in (a, b), with the boundary conditions M(a) = 0, M{b) = Sq. 

3 The case of an even number of given moments: 
solvability and a description of solutions. 

Consider the matrix moment problem ([1]) with I = 2d -\- 1, d € Z_|_. 
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Theorem 3.1 Let a matrix moment problem HP with I = 2d + 1, d G Z_| 
be given. The moment problem has a solution if and only if 



Td > 0; Td > 0, 
and there exist matrix solutions X, Y of matrix equations 



(73) 



( Sd+i \ 

Sd+2 
\ S2d+1 I 



VdY 



( -abSd + (a + b)Sd+i - Sd+2 \ 
-abSd+i + (a + b)Sd+2 - Sd+s 

\ -abS2d-i + (a + b)S2d - S2d+i ) 



(74) 



and for these solutions X, Y the following relation is true: 
X*TdX < -abS2d + (a + 6)52^+1 - Y*fdY. 



(75) 



Proof. Consider a matrix moment problem ([T|) with I = 2d + I, d TL+. 
It has a solution if and only if the moment problem with an odd number of 
moments 



x'^dMix) = Sn, n = 0, 1, 2d + 1; 



x^'^+'^dM{x) = S- 



(76) 



with some complex (A^ x N) matrix S2d+2 bas a solution. By ([5]) the solv- 
ability of the moment problem (j76p is equivalent to the matrix inequalities 



Td+i > 0, Td+i > 0. 



(77) 



If we apply to the latter inequalities the well-known block-matrix lemma 
(e.g. [171 P- 223]), we obtain that solvability of ([76]) is equivalent to the 
condition ([TH]) , existence of solutions X, Y of ([71]) and inequalities 

S2d+2 > X*TdX, S2d+2 < -abS2d + {a + b)S2d+i - Y*fdY. (78) 
Consequently, we get that the statement of the Theorem is true. □ 

Theorem 3.2 Let a matrix moment problem HP with I = 2d+l, d € Z+, be 
given. Suppose that conditions 1^7 3\ ), |7^[ ) and ( |75[ ) are true. For an arbitrary 
S2d+2 £ [X*TdX < —abS2d + {0'-\-b)S2d+i —Y*'^dY]! the set of all solutions 
of the matrix moment problem with an odd number of moments |76p we 
denote by V{S2d+2)- The set V of all solutions of the moment problem ([ip 
is given by the formula 



v = u, 



MS). 



(79) 



'Se[X*raX-abS2d + {'^+b)S2d+i-Y*TaY] 

The sets V{S) in ( [79[ ) for different S do not intersect. Each set V{S) in fyP] 
is parameterized by virtue of Theorem 
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Proof. The proof of (j79p follows obviously from the above considera- 
tions in the proof of Theorem 13.11 The sets V{S) in (|79p for different S do 
not intersect since the solutions in different sets have different {2d + 2)-th 
moments. □ 
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Truncated matricial moment problems on a finite interval: the 

operator approach. 

S.M. Zagorodnyuk 

In this paper we obtain a description of all solutions of truncated matricial 
moment problems on a finite interval in a general case (no conditions be- 
sides solvability are assumed). We use the basic results of M.G. Krein and 
I.E. Ovcharenko about generalized sc-resolvents of Hermitian contractions. 
MSG 2000: 44A60, 30E05. 
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